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Abstract 

An eight-potential-well order-disorder ferroelectric model was presented and 
the phase transition was studied under the mean-field approximation. It was 
shown that the two-body interactions are able to account for the first-order 
and the second order phase transitions. With increasing the random fields in 
the system, a first-order phase transition is transformed into a second-order 
phase transition, and furthermore, a second-order phase transition is inhib- 
ited. However, proper random fields can promote the spontaneous appearance 
of a first-order phase transition by increasing the overcooled temperature. The 

connections of the model with relaxors were discussed. 
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Relaxor ferroelectrics (relaxors) experience no macroscopic phase transition at zero elec- 
tric field.S Nevertheless, a ferroelectric phase can be induced by an external electric field at 
low temperatures. For example, in Pb(Mgi/3Nb2/3)03 (PMN), a first-order induced phase 
transition, from the pseudo cubic phase to a macroscopically polar phase, has been disclosed 
by Ye et a/..i A field-temperature phase diagram was proposed for PMN and a minimum 
threshold field strength Eth to induce the ferroelectric phase was determined.li The random 
interactions between polar microregions and the random electric fields existing in the sys- 
tem were usually adopted to interpret the special characteristics of relaxors. Glinchuk et 
al. showed that strong random fields can make the long-range order parameter disappear, 
and thus qualitatively explained the absence of the zero-field phase transition in relaxors.i 
However, they only considered a second-order phase transition and no external electric field 
was involved. It appears, therefore, necessary to investigate the effects of random fields and 
external fields on the first-order phase transition. 

A pseudospin mo del0 is usually used to interpret the phase-transition phenomena of 
order-disorder ferroelectrics in which the active ions have two equilibrium sites. In order 
to explain the first-order phase transition, four-spin coupling terms should be introduced 
and the four-body coupling strength is required to be larger than the two-body one.i On 
the other hand, PMN presents a disordered pseudo cubic structure with important atomic 
shifts around the special positions of the ideal perovskite structure,! i.e., the polar ions in 
PMN have more than two equilibrium positions. Therefore, a multi-direction "spin" system 
should be studied in order to understand the special properties of PMN. 

For the reasons mentioned above, we present an eight-potential-well ferroelectric model 
and investigate the the phase transition under various condictions in this paper. 

For PMN, the electric field-induced macropolar phase has trigonal symmetry with point 
group 3m,i which suggests that for each polar ion there are eight potential wells along {111}- 
equivalent directions. When an external field is applied along one of eight {lll}-equivalent 
directions, the Hamiltonian matrix of the order- disorder system with eight potential wells is 
assumed to be: 
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where E is the electric field strength, and po is the magnitude of the dipole when an ion 
locates in a certain well. A factor 1/3 is introduced in terms of Eq. (1) when the directions 
of the well are not parallel to that of the external field, fl is the tunneling frequency. Only 
the tunneling between wells with opposite directions are considered here. In ferroelectric 
systems, polar ions interact with each other and with different kinds of defects in the system. 
Under a mean field approximation, the interactions upon a certain ion may be represented 
by an equivalent field, i.e.. 



Epo — J h EextPO + ErandPO: 
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where (p) is the thermal average value of dipole moments, and J is the coupling energy. 
Only two-body interactions are involved here because a linear relation is assumed between 
(p) and E. E^xt is the applied external electric field. Erand is the internal random fields in 
the system which come from the direct interactions of dipoles, fields created by point charge 
defects and atomic ordering, etc.. E^and is assumed to have a Gaussian distribution: 

1 



P{Erand) 



27ra, 



■ exp 



_/?2 

rand 



(3) 



where (7^ is the distribution width. By using Eqns. (l)-(3), the order parameter (p) can be 
determined for any temperature. 

At first, the phase transition temperatures arc calculated when there is no external field 
and random fields, which are shown in Fig. 1. It can be seen that the long-range order 



can appear only at Q < 2 J. If the tunneling movement does not exist (fi = 0), the critical 
temperature is = where /c^ is the Boltzmann constant. An interesting characteristic 
of Fig. 1 is that both the second-order phase transition and the first-order phase transition 
exist for the present model. It is quite different from the case of two-direction pseudospin 
mo deli where no first-order phase transition exists if the number of coupled spins is two. In 
Fig.l, there is a tricritical point between the second- and the first-order phase transitions, 
which is determined as Qfri = 0.56 J and Ttri = 0.23 J/ ks- For a first-order phase transition, 
the overcooled temperature T_ decreases rapidly with increasing Q and vanishes at i7 = | J. 
When I J < r2 < 2 J, the long-range order cannot spontaneously appear in a cooling process, 
while a polar state can be induced at low temperatures. This is qualitatively consistent with 
the behavior of PMN.i 

Solid lines in Fig. 2 depict the phase transition temperatures versus external field when 
Q = J and (Tg = 0. At small external fields, the phase transition is the first-order type. 
And at large external fields, a second-order phase transition occurs. A tricritical point also 
appears here. For a first-order phase transition, the induced ferroelectric phase disappeares 
at the overheated temperature in a heating process. It is seen from the figure that the 
overheated temperature increases with increasing the field strength, whereas the zero-field 
depoling of an induced state occurs always at the temperature T^o (see Fig. 2), indepen- 
dently of the initial poling field. The existence of T^o and the tricritical point agrees with the 
experiments in PMN.i0 The curve of the overcooled temperature can also be interpreted as 
the curve of the threshold field strength Eth to induce a ferroelectric state at different tem- 
peratures. Eth decreases with decreasing temperature in our calculation. In experiments,!! 
however, Eth decreases firstly and then increases with decreasing temperature. The discrep- 
ancies come from the dynamical process, which is not considered in the calculation. At low 
temperatures, the system is frozen into a glass-like state.iS The dynamical process may be 
so slow that the induced polar phase cannot come into being even if the polar phase has 
smaller free energy. As a result, Eth increases at low temperatures (T < T^o). 

In Fig. 2, when the external field is larger than the tricritical point, there is no "real" 
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phase transition in the sense that the spontaneous breakdown of symmetry does not occur. 
Nevertheless, the curves of the average polarization (p) in cooling processes (see Fig. 3) 
indicate that (p) increases rapidly in a small temperature range but changes slowly at other 
temperatures when the external field is a little larger than the tricritical value (curve c in 
Fig. 3). Thus, a second-order "phase transition" may be approximately defined at the 
temperature where (p) changes most quickly. When the external field is much larger than 
the tricritical value (curve d in Fig. 3), (p) varies rather smoothly, and the phase transition 
is difficult to determine. The features of the polarization curves in Fig. 3 are consistent 
with experiments.0 

To explore the effects of random fields, the phase transition temperatures for the random 
fields (Te = 0.2J/po are also depicted in Fig.2 (dashed lines). Compared with the curves with 
no random field (solid lines), the critical temperature and the overheated temperature de- 
crease, and the range of the first-order phase transition gets smaller, i.e., the phase transition 
is compressed. It should be noted that the overcooled temperature increases when random 
fields exist. Thus the phase transition is promoted by the random fields in the sense that a 
spontaneous phase transition is easier to occur, which is in conflict with the common sense. 

The continuous effects of random fields on the phase transition temperatures are shown in 
Fig. 4 for Q = J. By enhancing random fields, the interval between the overheated and the 
overcooled temperatures shrinks gradually, and vanishes at a tricritical point. Beyond the 
tricritical point, the random fields make the critical temperature decrease, and finally inhibit 
the phase transition completely. The promotion effect of random fields on the spontaneous 
first-order phase transition is obvious in Fig. 4: the critical temperature keeps approximately 
constant while the overcooled temperature increases extraordinarily. This unusual effect 
seems to challenge the previous idea that the phase transition in PMN is inhibited by strong 
random fields. Perhaps the actual situation is that the random fields are not strong enough 
to promote a spontaneous phase transition. 

At last, a phase diagram on the Q — ae plane is given in Fig. 5. The Q—Ce plane is divided 
into three regions where a first-order phase transition, a second-order phase transition and 
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no phase transition can occur, respectively. The boundary hne between the first-order and 
the second-order phase transition regions is a tricritical hne that is composed of tricritical 
points. 

In summary, the phase transition of an eight-potential-well order-disorder ferroelectric 
system was investigated in this paper. The effects of the couplings, the external fields, and 
the internal random fields are discussed. The connections of the model with PMN were 
discussed. 

This work was supported by the Chinese National Science Foundation and State Key 
Program of Basic Research Development, and the authors would like to thank Dr. Jian-She 
Liu for helpful discussions. 
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FIGURES 

FIG. 1. Phase transition temperatures as functions of tunneling frequency with no random 

fields. The sohd line, the dashed line, and the dotted line represent the critical, the overheated, 
and the overcooled temperatures, respectively. The temperatures and the tunneling frequency are 
measured in units of J/ks and J respectively. 

FIG. 2. Phase transition temperatures as functions of external field when = J. The lines 
above tricritical points represent critical temperatures. Below tricritical points, lines on the left 
represent overcooled temperatures, and that on the right represent overheated temperatures. The 
solid and dashed lines correspond to the case of (Je = and 0.2J/po, respectively. External fields 
are measured in unit of J/po and temperatures in unit of J/ks- 

FIG. 3. Average polarization {p) (in unit of po) as function of temperature (in unit of J/ks) 
in cooling processes at = J and = 0. Curves a, b, c, and d correspond to the external field 
E^^t = 0.075, 0.095, 0.115, and 0.15 J/pq, respectively. 

FIG. 4. Phase transition temperatures (in unit of J/ks) versus random fields width (in unit 
of J/pq) when J7 = J. The solid, dashed and dotted lines represent the critical, overheated, and 
overcooled temperatures, respectively. 

FIG. 5. Phase diagram on the — cxe plane, a, /3, and 7 designate the regions where the 
first-order, the second-order, and no phase transition can occur. f2 and Ue are measured in units 
of J and J/pO) respectively. 
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